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Third Zemach Moment of the Proton 
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Modern electron scattering experiments have determined the proton electric form factor, Ge p (Q 2 ), to 
high precision. We utilize this data, represented by the different empirical form factor parametrizations, 
to compute the third Zemach moment of the proton charge distribution. We find that existing data rule out 
a value of the third Zemach moment large enough to explain the current puzzle with the proton charge 
radius, determined from the Lamb shift in muonic hydrogen. This is in contrast with the recent paper 
of De Rujula. We also demonstrate that the size of the third Zemach moment is largely governed by the 
fourth moment of the conventional charge distributions, (r 4 ), which enables us to obtain a rigorous upper 
bound on the magnitude of the proton's third Zemach moment. 
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Pohl et al. recently reported a very precise measurement 
of the Lamb shift in muonic hydrogen, that leads to a pro- 
ton radius of [1] 



2\l/2 



0.84184 ± 0.00067 fm, 



(1) 



which is five standard deviations away from the CODATA 
compilation of J2] 

( r 2)V2 = 0.8768 ± 0.0069 fm. (2) 

There has been much speculation about the origin and pos- 
sible consequences of the puzzling difference between the 
results of muonic and electronic hydrogen. 

Even more recently De Rujula [3] pointed out that the 
discrepancy between the two results can be removed if the 
proton's third Zemach moment is very large. He finds that 
a value of 



(2) 



36.6±6.9fm J 



(3) 



combines the impressive experimental results from the 
CODATA compilation and the recent muonic Lamb shift 
data in a consistent manner. This result is more than 13 
times larger than the experimental extraction of Friar and 
Sick 0j, who use electron-proton scattering data to deter- 
mine 
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2.71 ±0.13 fm 3 



(4) 



It is noteworthy that the length scale associated with Eq. (O 
is about four times the proton root-mean-square radius. 
Such a large value is supported in Ref. [3] by a "toy model" 
of the proton electric form factor, namely 
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where s 2 = sin 6, c 2 = cos 2 9 and D = c 2 M 2 + s 2 
The CODATA charge radius result and the large Zemach 
moment of Eq. d3j> can be accommodated for a range of 
mass parameters and mixing angle. For s 2 = 0.3 the mass 
scales are typically m ~ 15MeV and M ~ 750 MeV [0]. 
For these parameters we have plotted Eq. (01 as the dotted 
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FIG. 1 : The proton G e p form factor from the parameterizations 
of Kelly [31 (solid line), fit II from Alberico et al. |6] (dashed 
line) and the dipole form factor of Eq. © is given as the dot- 
dashed line. The "toy model" of Ref. [J], with sin 2 6 = 0.3, 
m = 15 MeV and M = 750 MeV, is illustrated as the dotted line. 



line in Fig.Q] where it is compared with a dipole form fac- 
tor and the empirical parameterizations from Refs. 15|,|6|]. 
The severe disagreement between the form factor of Eq. (f5]l 
and the realistic form factors of Refs. @] demonstrate 
that this "toy model" is not a viable representation of the 
data. 
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parameterization 


(^>(2) [fm 3 ] 


(r 2 p ) [fm 2 ] 


(r 2 ) 1/2 [fm] 


(4) [fm 4 ] 


(r 3 p )f% [fm 3 ] 


Dipole (Eq. ©) 


2.023 


0.658 


0.811 


1.083 


3.029 


Kelly [5J 


2.524 


0.744 


0.863 


1.619 


4.067 


Alberico et al. [6] 


2.545 


0.750 


0.866 


1.623 


4.097 



TABLE I: Values of the third Zemach moment (rp}^ 2 ), the mean-square charge radius (r 2 ,), the charge radius {r" 2 ,) 1 / 2 and the fourth 
moment of the charge distribution (r^), for each of the three proton form factor parameterizations considered. The last column 
containing (rp) 1 ^) represents the upper bound on the third Zemach moment of the proton given by Eq. dl5h with /i = 1 GeV. 



Nevertheless, it is reasonable to ask if the use of a realis- 
tic extraction of Ge p leads to a value of the Zemach radius 
very different from that obtained from the typically used 
dipole form, namely 



G E pole (Q 2 ) 
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where A 2 = 0.71 GeV 2 . This form has historical valid- 
ity in describing early data and was explicitly assumed in 
the recent experimental analysis of Pohl et al. [1]. How- 
ever, several experiments (see for example the review of 
Ref. 170) have found that G Ep actually falls faster than 
the dipole. Therefore, we use two different more recent 
parametrizations from Refs. Jli@]. These parametrizations 
take the general form 



1 + Ett^T" 
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Q 2 

4 m 2 : 
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with proton mass labeled by m p and the various coeffi- 
cients are given in Refs. JUta]. Note, we use fit II for the 
Alberico et al. [6] parameterization. These empirical re- 
sults for the proton G Ep form factor are illustrated in Fig.Q] 
The third Zemach moment is defined by 



where 



(rp( 2 ) = J d i rr s p 2 {r), 

p 2 (r) = j d 3 r'p p (r')p p (\r'-r\), 
-J (2vr)3 6 G * (qj - 
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The conventional proton charge density, p p (r), is defined 
by 

P P (r) = J 703^^ G Ep (q 2 ), (10) 

with the charge radius given by 

(r 2 ) = J d 3 rr 2 p p {r). (11) 

The results of our numerical evaluations are presented 
in Table U Observe that the values of (?"p}(2) ar e smaller 



than that of Eq. (01 by approximately a factor of 15. How- 
ever, the agreement with the empirical result of Ref. 0] 
(see Eq. ©) is fairly good. The G Ep parameterizations of 
Refs. yL |6|] fall faster with increasing Q 2 than the dipole 
form factor (see Fig. [Q, leading to coordinate space dis- 
tributions of greater extent. Therefore, the values of the 
third Zemach moment (fp)(2) and the mean-square radius 
(r 2 ) both increase relative to the values obtained from the 
dipole form factor. 

Examination of the integrand in Eq. dSJ for each of the 
three form factor models, reveals that 80% of the strength 
for the Zemach moment comes from the domain where Q 2 
is less than 1 GeV 2 and 50% from Q 2 values less than 
0.25 GeV 2 . Therefore, the quoted values of the Zemach 
moment given in Table [I] are insensitive to the model de- 
pendent extrapolation to infinite Q 2 . In the experimental 
extraction of the charge radius by Pohl et al. QJJ], the third 
Zemach moment is suppressed relative to the charge radius 
by a factor of ~ 570, therefore the variation of the Zemach 
moment observed in Table U has negligible impact on the 
extracted value of the charge radius. 

It is possible to obtain a rigorous upper bound on the 
third Zemach moment of the proton using the result [0] 
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G Ep (q ) + \( r p) 



(12) 

The integrand in Eq. (fT2l is a monotonically decreasing 
positive definite function on the integration domain, with a 
maximum at q = given by 
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For moderate to large q the integrand in Eq. (fT2l is com- 
pletely saturated by the last two terms, that is 
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Empirically we know that at q = lGeV the above ex- 
pression represents ~99.5% of the integrand in Eq. (fT2l . 
This is illustrated in Fig. |2] for the dipole form factor 
of Eq. (0, where we plot the integrand of Eq. (fT2l and 
Eq. (fT4l . The integrand of Eq. dT2l) coincides with Eq. (fT4l 
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FIG. 2: Plot of the integrand in Eq. d 12b (solid line) and the dom- 
inate piece at moderate to large q given by Eq. d 14-b (dashed line). 

f° r Q ^ 1 GeV, this is also the case for the empirical form 
factors of Refs. 1110] and the model of De Rujula. 

Therefore, the third Zemach moment has a rigorous up- 
per bound given by 




where the first term is the area of the rectangle defined by 
the integrand upper bound of Eq. (fT3l l and ^ q ^ \i GeV, 
whereas the second term results from integrating Eq. ([T4il 
from q = //GeV to infinity. For the reasons discussed 
above we expect that fj, = 1 GeV is a conservative choice 
for the evaluation of the upper bound given in Eq. ( fl~5l) . It 
is therefore apparent that if a model reproduces the proton 
charge radius, a large third Zemach moment can only occur 
if (rp) is also large. Typical empirical values are given in 
the second last column of Table[]] whereas the model of De 
Rujula, expressed in Eq. (O, gives (rf) = 1849 fm 4 . The 



upper bound for the third Zemach moment of the proton 
obtained from Eq. (fl~5l . with fx = lGeV, is given in the 
last column of Table U where it is labeled by (r^)™". 

The net result is that the published parametrizations, 
which take into account a wide variety of electron scatter- 
ing data, cannot account for the value of the third Zemach 
moment found in Ref. yfl. It is likely that further investiga- 
tions of low Q 2 data will lead to improved parametrizations 
for Ge p S- However, enhancing the Zemach moment sig- 
nificantly above the dipole result is extremely unlikely. We 
believe that the resolution of the current puzzle regarding 
the proton radius, determined from the muon and electron 
Lamb shifts in hydrogen, lies in a direction other than that 
suggested by De Rujula. 
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